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THE  RESISTANCE  OF  A  FLAT  PLATE  LOCATED  NORMAL 
TO  THE  FLOW  OF  GREATLY  RAREFIED  GAS 


V.  A.  Perepukhov 
(Moscow) 

Given  in  the  work  is  a  method  of  calculation  of 
aerodynamic  properties  of  a  plate  of  arbitrary  form  with 
an  arbitrary  value  of  the  coefficient  of  accommodation  in 
the  flow  of  greatly  rarefied  gas,  taking  into  account 
first  collisions  between  molecules  of  reflected  and 
incident  flows.  As  an  example  calculg tions  of  the 
aerodynamic  properties  of  a  square  are  given. 

1.  Introduction.  In  works  [1-8]  there  was  x  /estigated  the 
streamline  flow  of  bodies  of  different  form  possessing  a  coefficient 
of  accommodation  =  0  in  the  range  of  the  so-called  hyperthermal 
theory,  taking  into  account  first  intermolecular  collisions.  Let  us 
recall  the  basic  positions  of  this  theory. 

At  very  large  macroscopic  speeds  of  flows  of  greatly  rarefied 
gas,  the  characteristic  speed  in  the  flow  is  the  mean  thermal  velocity 
of  the  molecules  VT  —  V 2JiTm  and  n'ot  the  speed  of  sound,  as  it  was 
in  the  hypersonic  theory.  Therefore,  instead  cf  the  dimensionless 
number  M  the  new  number  Sm=  UJY 2RT„  is  introduced. 

Let  us  assume  that  molecules  of  the  flow  incident  on  the  body 


possess  macroscopic  velocity  Uro,  much  greater  than  the  mean  thermal 
velocity,  i.e.,  o„  I .  Not  introducing  gross  errors  into  further 


calculations,  it  is  possible  to  assume  that  the  incident  flow  consistr 

of  molecules  moving  in  parallel  to  each  other  at  macroscopic  speed 

Um  (such  an  assumption  is  acceptable  when  £  5). 

Let  us  assume  that  the  temperature  of  the  body  is  low  because 

of  heat  removal  so  that  T  ~  Tro;  since  the  coefficient  of  accommodation 
* 

ctA  ~  0  the  speed  of  the  reflected  molecules  has  an  order  of  thermal 
velocity  of  molecules  of  incident  flow,  which  permits  disregarding 
the  latter  at  the  time  of  collision  of  the  molecule  of  incident 
flow  with  the  reflected  speed  as  compared  to  the  speed  U  . 

In  works  carried  out  earlier  there  is  proposed  the  following 
method  of  calculation  of  aerodynamic  properties  of  a  body.  On  the 
surface  of  the  body  an  element  of  the  surface  dF  was  selected,  and 
there  was  calculated  a  change  of  some  molecular  criterion  arriving 
on  the  element  dF  because  of  collisions  of  molecules  of  incident 
flow  with  molecules  reflected  from  the  whole  surface  of  the  body 


rw  =  — ri(.)/dF, 

where  —  flow  of  the  j-th  molecular  sign  on  element  dF,  H0^.  — 

flow  of  the  j-th  molecular  criterion  on  element  dF  during  freely 
molecular  streamline  flow,  —  flow  of  the  j-th  molecular  sign 

arriving  on  the  element  dF  due  to  collision  of  molecules  of 
incident  flow  with  the  reflected,  -  flow  of  j-tb  molecular  sign 

which  because  of  collisions  of  molecules  of  incident  flow  with  the 
reflected  does  not  fall  on  the  element  dF. 

In  this  work  the  solution  of  another  problem  with  the  help  of 
a  calculation  different  from  the  preceding  method  is  proposed. 

2.  Formulation  of  problem.  Let  us  consider  the  streamline  flow 
of  a  flat  plate  of  arbitrary  form  in  the  plan  possessing  an  arbitrary 
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coefficient  of  accommodation  and  arbitrary  temperature  Tw  and 
a  flow  of  greatly  rarefied  gas. 

The  lav/  of  reflection  is  best  to  select  on  the  basis  of  the 
experiment,  but  since  up  to  now  there  has  been  no  reliable 
experimental  data  then  we  will  select  the  classical  law  of 
reflection,  diffuse  the  Maxwellian  function  of  distribution  of 
reflected  molecules  and  temperature  TQ  (i.e.,  with  energy  £Q) : 


/o  =  ho  (~VVh*‘\ 


ho  = 


m 

UTc  ‘ 


The  presence  of  the  arbitrary  coefficient  of  accommodation 
leads  to  the  fact  that  at  the  moment  of  collision  of  the  molecule 
of  incident  flow  with  the  reflected  speed  the  latter  cannot  be 
neglected,  which  in  its  turn  leads  to  the  complication  of  the 
integral  of  collisions.  Thus  as  earlier,  v/e  will  assume  that  the 
molecules  are  solid  balls  with  a  diameter  a. 

Of  all  the  forms  of  possible  collisions  of  molecules  among 
themselves  belonging  to  different  flows,  v/e  will  consider  only  the 
collisions,  of  molecules  of  incident  flov/  with  reflected  flow, 
considering  the  effect  of  damping  of  the  x’eflected  flow  on  the 
incident. 

Taking  into  account  those  difficulties  which  are  found  in  the 
solution  of  the  preceding  problems  and  the  laboriousness  of  the 
method  of  calculation,  in  this  problems  v/e  'will  use  a  somewhat 
different  method  founded  on  the  calculation  of  the  influence 
function  for  some  molecular  criterion  of  one  element  of  the  surface 
on  another  [6]. 

Let  us  examine  the  small,  as  compared  to  the  surface  of  the 
body,  element  of  the  surface  dF^  and  calculate  what  is  influence  of 
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molecules  reflected  from  this  element  on  aerodynamic  properties  of 

another  element  of  the  surface  dF,.  By  knowing  the  influence  function 

J 

of  element  dFi  on  the  remaining  elements  of  the  plane,  one  can 
determine  (in  virtue  of  linearity  of  the  problem)  the  influence  of 
all  elements  entering  into  the  plate  of  a  given  configuration  on  the 
element  of  plate  dF^,  and  then  by  adding  the  flow  of  the  corresponding 
molecular  sign  on  the  whole  plate  we  will  obtain  the  total  flow  of 
some  molecular  sign  for  the  whole  surface  of  the  body. 

5.  Diagram  of  solution  of  the  problem  on  a  computer.  We  will 
solve  this  problem  with  the  application  of  the  Monte-Carlo  method, 
v.'hich  in  preceding  works  was  used  for  the  calculation  of  the  integral 
of  collision. 

Let  us  fix  in  space  near  the  element  dFi  an  arbitrary  point 
A(X,  Y,  Z),  then  select  an  velocity  vector  arbitrary  in  magnitude 
of  reflected  molecules  vQ;  its  direction  is  determined  by  coordinates 
of  the  element  and  \oint  A(X,  Y,  Z) . 

Let  us  assign  further  the  impact  parameters  ¥  and  e  (Fig.  1) 
with  respect  to  the  direction  of  speed  g  =  Ura  -  vQ  and  determine 
by  it  the  directing  cosines  of  the  line  of  centers  at  the  time  of 
collision  of  the  reflected  molecule  with  the  molecule  of  incident 
flow  (in  the  system  of  coordinates  connected  with  the  element  dFi). 
Then  the  quantity  of  collisions  in  the  element  of  volume  at  point  A 
is  such  that  the  reflected  molecules  possess  speed  in  the  element 
dV,  and  the  line  of  centers,  found  in  the  solid  angle  dT,  will  be 
equal  to 
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Fig.  1. 


where  g1  —  Vo+  UJo  +  2u0Um.cos  [i 

X  =  r  sin  n  cos  p,  Y  =-r  sin  n  sin  P ,  Z  —  r  cos  p. 

The  expression  for  nQ  can  be  obtained  from  the  law  of  preservation 
of  the  number  particles  on  the  surface  of  the  plate 

n0  =  2nmU oo  (nAo)^*. 

The  direction  cosines  of  the  line  of  centers  with  respect  to 
the  system  of  coordinates  connected  with  element  dF^  can  be  \*ritten 
thus 

l  (a:)  =  cos  ¥  cos  P  sin  a  4-  sin  ¥  cos  a  cos  P  cos  e  — sin-¥  sin  e  sin  p; 
m  (y)  =  cos  'F  sin  p  sin  a  +  sin  ¥  cos  a  sin  p  cos  e  +  sin  ¥  sin  e  C03  P; 
n  (z)  =  cos  ¥  cos  a  —  sin  a  sin  ¥  cos  e, 


where 


The  variables  of  problem  have  the  following  range  of  variation: 


0;<  V  <  a/2,  0  <  I*  <  n/2,  0  <  e  <  2a, 
0  <  P  <  2n,  0  <  r  <  oof  0  <  v0  <  co. 


In  the  expression  for  the  number  of  collisions  the  recording 

_p 

used  for  the  solid  angle  r  dF  cos  u  is  accurate,  strictly  speaking, 
when  V'dF  <r  ,  i.e.,  with  a  tendency  r  —  0  dF  it  should  approach 
zero  faster.  However,  in  our  case  dF  is  the  fixed  magnitude  and 
with  a  tendency  r  -*•  0  the  solid  angle  is  equal  to  co.  in  order  that 
this  not  be  true,  it  is  possible,  not  exceeding  the  bounds  of 
accuracy  of  the  method,  to  record  the  solid  angle  in  the  following 


form: 


■  /  (r)  cos  p,  where 


/(r>- 

^  (+<»•) 


All  the  variables  of  the  problem  for  every  collision  can  be 

selected  equiprobable  within  limits  of  a  change  of  each  variable. 

Above  it  v/as  stated  that  for  every  element  of  surface  dF .  it  is 

3 

necessary  to  determine  two  forms  of  flow  of  the  n-th  molecular  sign 

with  the  index  (+)  and  (-).  Flows  with  the  (+)  index  correspond  to 

flows  arriving  on  the  element  dF .  due  to  collision  of  the  reflected 

J 

and  incident  molecules  at  an  arbitrarily  selected  point  of  space 

(V  Ym,  Zm).  Flows  with  the  (-)  index  correspond  to  those  flows 

of  molecular  signs  which  do  not  attain  the  element  dF .  due  to 

J 

collisions  of  reflected  and  incident  molecules  at  points  (X  .  Y  ,  Z  ), 

v  m*  m  nr 

i,e.>  at  points  located  above  the  element  dF.. 

J 

If  we  want  to  determine  the  flow  of  some  molecular  sign 

arriving  on  element  dF,.  as  a  result  of  dispersion  of  molecules  of 

J 

incident  flow  on  reflected  flow  from  the  element  dF.,  we  should 

J 

multiply  the  number  of  collisions  by  the  appropriate  weight 
Bn(±)*  Por  flow  mas s ,  the  flow  of  normal  pulse,  and  the 
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'  WB'  ' 


s  |  vi  '-.'■"TVY-’-  ’ ' 


energy  flow  expressions  for  and  ^n(-)  ^ave  following  form: 

/?1(+)  =  m,  ll2U)  =  tnv[l,i\  n3M  =*.„*{ i.i)( 

^?K-)  =  m,  =•  mUc o,  =  -j-Jl,, 

where  the  indices  1 ,  2,  3  below  correspond  co  the  above-mentioned 
flows,  and  indices  (1,  2)  above  correspond  to  two  molecules 
appearing  as  a  result  of  the  collision.  It  is  possible  to  obtain  an 
expression  for  speeds  of  molecules  after  collision  depending  upon  the 

value  of  these  speeds  up  to  ’’.h.*  •'ollision  and  impact  parameters 

»»>  =  v0  sin  fx  cos  p  —  la,  yd>  —  *  sin  {i  sin  P  —  tna,  yd>  =  v0  cos  n  —  no,. 

v{*:  =  la,  yd>  =  ma,  yd>  =  —  U„,  -4-  na, 

to  —  /-v  sin  u  cos  P  +  niv 0  sju  ft  sin  p  +  n  ( Um  +  v0  cos  fi), 

i ,  m,  n  are  direction  cosines  o^  the  line  of  centers  at  the  time  of 
collision  in  the  system  of  coordinates  connected  with  the  element  of 


After  the  collision  of  the  trajectory  of  the  molecules  one  can 

determine  the  direction  cosines  in  the  system  of  coordinates  with 

the  center  at  point  A(X  ,  Y  ,  Z  ) : 

'  m  m  m7 


)<« 

*"X 

KM 

v(,) 

Vy 

Vv 

v<«  _ 

vv 

]/v;<»  +  *(.)  +  vf(«) 

V 

t(.) 

•  T?>=- 

yd) 

V1 

— - - -  j 

vz 

r— - : - : —  t 

V vl{l)  *r  V0  +  vi(,) 


r(2) 


(1) 


if  v;  ’  <  0  or  v'  '  <0,  then 

u  z 

id) 

X1  =  Xm  +  ^TZm, 
■  • 

V<1) 

Y  l  —  Fm  H - Zm, 


x(” 

X*  =  Xm  +  -± -zm, 
Y,  =  Ym  +  -fcZm, 

T t 


i.e.,  after  the  given  collision  the  molecule  will  arrive  in  elements 
dF^  and  dFg  with  the  corresponding  speeds  v^1)  and  v^2). 

Thus  it  is  possible  to  construct  completely  the  form  of  the 
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influence  function  with  weight  Rn  °f  element  dP^  on  any  element  of 
plane  dF.,  after  which  the  calculation  of  integral  aerodynamic 

d 

properties  is  not  difficult. 

4.  Results  of  calculation.  A  calculation  was  produced  for  the 
influence  functions  corresponding  to  the  particle  flux,  bo  the  flux 
of  the  normal  pulse,  and  to  the  energy  flux.  With  this  the  influence 
was  determined  of  the  square  element  on  a  certain  quantity  of  other 
elements,  so  that  as  a  result  of  the  calculation  it  became  to 
determine  the  aerodynamic  properties  of  a  flat  plate  arbitrary  form 
inscribed  into  square  (il  x  11) . 

A  calculation  was  produced  for  the  following  values  of  the 
coefficient  a*  =  1,  1/2,  1/52,  1/128,  where 

a*A  =  8RTJUI- 

The  corresponding  influence  functions  with  a  different  value  of 
are  given  in  Tables  1-3.  Results  of  the  calculation  when 

■Xr 

aA  =  1/128  practically  do  not  differ  from  results  of  the  calculation 
when  »  1/32.  In  order  to  obtain  a  correction  for  some  flow  of 
molecular  sign  arriving  on  element  dP^  it  is  necessary  to  add  the 
influence  functions  for  the  corresponding  criteria 

ACD  jkJ  =  dF{  ^  W'ii1  dFj? 

J 

o 

Here  the  Knudsen  number  Kn^  =  (2Tra‘_n00)  ,  since  the 

characteristic  dimension  is  equal  to  1  (unit  of  length  is  the  side 
of  the  unit  square  element) . 


Wjf- 


which  corresponds  to  the  fluxes  of  molecules,  normal  pulse,  and 


energy. 
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The  flow  of  the  k-th  molecular  criterion  arriving  on  elemei  t  aF^ 


can  be  written  in  the  following  form: 


(D{x,=  A(DiM, 

where  <!>(,«  —  is  the  value  of  the  corresponding  flov;  of  molecular 
criterion  in  freely  molecular  conditions  of  streamline  flov;.  For 
the  flov;  of  normal  pulse  and  energy  flow  it  is 


-Si- 


necessary  to  consider 


Fig.  4. 


the  reactive  value  of  the  corresponding  flows. 

As  an  example  theie  is  produced  a  calculation  of  the  aerodynamic 
properties  of  a  square  plate  with  side  d.  For  an  appraisal  of  the 
accuracy  of  the  method  a  comparison  is  made  with  the  calculation 
produced  in  [8]  when  -+  <x>  (Sq  co) . 

Figure  2  gives  curves  (s^lid  line  corresponds  to  the  distribution 
along  the  diagonal,  dashed  line  —  along  the  side)  for  correction  to 
the  flow  particles 

Ni  “  A'<  Kt\aJnmUaS0. 

Figures  3-4  give  curves  for  corrections  to  the  flow  of  the 
normal  pulse  and  to  the  energy  flow  without  taking  into  account  the 
reactive  flow 


and  for  the  energy  flow 

5T=Tt[1-i+^:(i#'7*)+w(^),]»  ' 

the  corresponding  curves  are  given  in  Pigs.  5  and  6  where 


LP  ^N,X?-+  P{  and  Lb  =  A  Af<- Et. 


For  the  total  fluxes  of  the  mass,  normal  pulse,  and  energy 
the  corresponding  expressions,  taking  into  account  the  jet  stream 

O 

(disregarding  terms  (S0/Knco)  and  of  higher 
order)  are  obtained 

I 

(the  total  flow  of  the  mass  is  referred  to 

»w> » 


Pig.  7. 


Cx  —  2  + 


Y~ 

So  ( 

So  ~ 

"Knro  [ 

(?*£+?*) 


(the  total  flow  of  the  normal  pulse  is  referred  to  1/2  mn^U^F), 

(the  total  flov;  is  referred  to  1/2  mn^U^F) .  Figure  7  gives  curves 
plotted  depending  upon  Sq  for 

{  .  ,1  .  t  *0  t  l 


Ur 


==■ .  Kn^  =  {2iX!5s/i00rf]_a. 


V2M 

In  conclusion  the  author  expresses  sincere  gratitude  to  M.  N. 
Kogan  for  a  discussion  of  the  results  and  interest  toward  the  work. 


Submitted 
26  July  1964 


FTD-NT-67-25 


-15- 


LITERATURE 


1.  Hyneman,  M.  K.  11 K  teorii  lobovogo  soprotivleniya  v  sil'no 
razrezhennykh  gazakh"  (The  theory  of  drag  in  greatly  rarefied  gases) 
Mekhanika,  Sb.  perev.  i  oboz.  in.  period,  lit.,  No.  3,  1951. 

2.  Lunts  and  Lyubovskiy.  "Obtekaniye  prepyatstviya  svobodno 
molekulyarnym  potokom"  (Streamline  flow  of  an  obstacle  by  freely 
molecular  flow)  Mekhanika,  Sb.  perev.  i  oboz.  in.  period,  lit..  No.  2, 
1958. 


5.  Lyu.  "0  soprotivlenii  ploskoy  plastinki,  obtekayemoy 
pod  nulevym  uglom"  (The  resistance  of  a  flat  plate  of  streamlined 
at  a  zero  angle)  Mekhanika,  Sb.  perev.  i  oboz.  in.  period,  lit.. 

No.  1,  I960. 

4.  Perepukhov,  V.  A.  11 0  soprotivlenii  ploskoy  plastiny  v 
potoke  sil'no' razrezhennogo  gaza"  (The  resistance  of  a  flat  plate 
in  a  flov;  of  greatly  rarefied  gas)  Zh.  vychisl.  matem.  i  matem. 
fiz.,  Vol.  1,  No.  4,  1961. 

5.  Kogan,  M.  N.  "0  giperzvukovykh  techeniyakh  razrezh  ennogo 

gaza"  (Hypersonic  flows  of  rarefied  gas)  Prikl.  matem.  i  me khan. , 
No.  3,  1962.  ' 

6.  Kogan,  M.  N.  "Teorema  obratimosti  dlya  techeniy,  blizkikh 
k  svobodnomolekulyarnym"  (Reciprocity  theorem  for  flows  close  to 
free-molecular)  Dokl.  AN  SSSR,  Vol.  144,  Ho.  6,  1962. 

7.  Perepukhov,  V.  A.  "Obtekaniye  ploskoy  plastiny,  raspolo- 
zhennoy  pod  nulevym  uglom  ataki,  potokom  sil'no  razrezhennogo  gaza" 
(Streamline  flow  of  a  flat  plate  located  at  a  zero  angle  of  attack 
by  a  flow  of  greatly  rarefied  gas).  Zh., vychisl.  matem.  i  matem. 
fiz.,  Vol.  No.  5,  1963. 

8.  Friedlander,  0.  G.  "0  soprotivlenii  ploskoy  plastiy, 
perpendikulyarnoy  giperzvukovomu  potoku  razrezhennogo  gaza" 

(The  resistance  of  a  flat  plate  of  perpendicular  to  the  hypersonic 
flow  of  rarefied  gas).  Zh.  prikl.  mekhan,  i  tekhn.  fiz.,  No.  3, 
1965. 


FTD-MT-67-33 


-14- 


